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ABSTRACT 

The HARPS spectrograph is showing an extreme stability close to the ms“^ level over more than 
ten years of data. However the radial velocities of some stars are contaminated by a spurious one-year 
signal with an amplitude that can be as high as a few ms“^. This signal is in opposition of phase 
with the revolution of Earth around the Sun, and can be explained by the deformation of spectral lines 
crossing block stitchings of the CCD when the spectrum of an observed star is alternatively blue- and 
red-shifted due to the motion of Earth around the Sun. This annual perturbation can be supress by 
either removing those affected spectral lines from the correlation mask used by the cross correlation 
technique to derive precise radial velocities, or by simply fitting a yearly sinusoid to the RV data. 
This is mandatory if we want to detect long-period low-amplitude signals in the HARPS radial 
velocities of quiet solar-type stars. 


Subject headings: instrumentation: detectors - instrumentation: spectrographs - methods: data analysis - techniques: 
radial velocities - planets and satellites: terrestrial planets 


1. Introduction 


HARPS is a fiber-fed, cross-dispersed echelle spec¬ 
trograph. The instrument is fed by two fibers, one for 
the target and the other for a reference lamp or sky. 
The spectrograph re-images the two fibers on a mosaic 
of two 2k4 CCDs, where two echelle spectra of 72 or¬ 
ders are formed. The resolution of the spectrograph is 
R= 115,000 with a resolution element oversampling of 


3.2 CCD pixels (for more information, see Mayor et al. 
20031 IPepe et al.|2002| i. 


To prevent any large drift in radial velocity (RV), 
HARPS has been designed to reach an extreme 
long-term stability of 0.01 mbar in pressure and 
0.01 K in temperature. The tiny RV instrumental 
drifts related to residual pressure and temperature 
changes are measured by a Thorium-Argon refer¬ 
ence lamp (recently upgraded to a Fabry-Perot in¬ 
terferometer calibration source, [Wildi et al.|2011| |. 


Due to those exceptional characteristics, HARPS, 
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and now its copy HARPS-North ^Cosentino et al. 


20121 , has been shown to be the most precise spec¬ 


trographs to search for extra-solar planets. As an ex 
ample, the HARPS RV of the magnetically quiet stars 
rCeti (HD10700) and HD85512 show a standard de¬ 
viation of 0.92 m s“^ and 1.05 m s“^ over more than six 
years of data, respectively (see Table 2 in |Pepe et ak 
( [2011 ) for more examples). 


The precision of the HARPS spectrograph might be 
even below the ms“^ level, as we know that RV mea¬ 
surements are contaminated by different kinds of stel¬ 
lar signals that induce variations that can be as high as 
a few m s“^. To our current knowledge, these stellar 
signals can be differentiated in four categories: 


• stellar oscillations produced by pressure waves 
propagating in the convective zone of the star 
([Dumusque et al.|2011at |Arentoft et al.||2008} 

[Kjeldsen et al.|2ofe I 


• stellar granulation that is the result of convection 
at the stellar surface ( [Dumusque et al.|2011a| ) 

• short-term stellar activity induced by stellar ro¬ 
tation in the presence of stellar active regions 


( Dumusque et al. 12014 Boisse et al.|2012t[M^ 
|unier et al.|2010t|Saar&Donahue|1997| > 
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• long-term stellar activity induced by stellar mag- 
netic cycles ([Meunier & Lagrange 2013[ Du- 


musque et al.|201 


In this paper, we present a new type of perturbing 
signal that is not induced by the stars themselves, but 
by the design of the CCD detector. In Section we 
highlight the presence of this signal in the RVs of sev¬ 
eral stars observed with HARPS. We then demonstrate 
the origin of this signal in Section and present a so¬ 
lution to correct for this signal in Section]^ Section 
concludes the paper. 


2. A one-year signal present in HARPS RV data 

Several stars that have been intensively observed 
with HARPS during several years show a low-amplitude 
signal at a period of one year. This is for example the 
case of HD128621 (aCenB), HD1461, HD154088, or 
HD31527. Note that other stars are also affected by 
this annual effect. 


2.1. HD128621 


In Dumusque et al. ( 2012 1, the authors present a de¬ 
tailed analysis of the RV of HD128621 and detect the 
presence of an Earth-mass planet orbiting at a short pe¬ 
riod of 3.24 days. In addition to this planetary signal, a 
detection of the binary RV drift induced by HD 128620 
(aCenA), a RV variation induced by active regions, 
and a long-term drift due to the stellar magnetic cycle 
are reported. As explained in |Lindegren & Dravlii^ 
(|2003|) and shown in [Dumusque et ~ ( |2011b ) and 


Meunier & Lagrang^ ( |2013 1, a stellar magnetic cycle 


induces a correlated effect in RV because active re¬ 
gions strongly suppress stellar convection. By fitting 
a 3rd order polynomial to the raw RVs between 2008 
and 2011 to remove simultaneously the binary and the 
magnetic cycle effect, a yearly signal is detected in the 
residuals (Figure top left panel). This one-year sig¬ 
nal cannot be detected in the Dumusque et al. ( 2012| ) 
published RV data, because this effect was discovered 
at the time, and was corrected by using the method de¬ 
scribed in Sec. lU 


The yearly signal is extracted from the data by fit¬ 
ting a sinusoid with a 365.25-day period, in addition to 
the 3rd order polynomial. This yearly signal that has 
an amplitude of 2.83 + 0.47 ms“^ is shown in the top 
panel of Figurej^ The bottom panel of the same figure 
shows the barycentric Earth radial velocity (BERV), 
i.e. the velocity of Earth in the direction of HD128621. 


The stellar velocity is the relativistic sum of the raw 
RV as measured by the spectrograph and the BERV. 
The BERV varies with a period of one year because 
Earth orbits the Sun, and its maximum amplitude de¬ 
pends on the position of the star in the sky, but can 
be as large as +30kms"^. This velocity variation can 
be translated to a shift of the spectrum on the HARPS 
CCD of +37 pixel^ Although the strongest signal 
in the BERV is the one with the one-year period, 
other signals with different periods are relevant for 
precise Doppler work at the nis“^ precision. Ex¬ 
amples of these other signals are the rotation of the 
Earth around its axis, the revolution around the 
Sun of the solar system planets or relativistic effects 
(see e.g. [Wright & Eastman|2014t for a detailed re¬ 
view). In this study, we focus on the yearly signal 
because it is the only one that seem to be detected 
in the data. 

The fact that an anti-correlation is observed be¬ 
tween the BERV for HD 128621 and the one-year sig¬ 
nal detected in the RVs (Figure top left panel) in¬ 
dicates that the two yearly signals are in opposition of 
phase. The BERV is the velocity of the Earth in the 
direction of the star, and the RV is the measure of the 
stellar velocity in the direction of Earth. Therefore this 
opposition in phase suggests that the one year signal 
detected in the RVs is induced by the orbit of Earth 
around the Sun and is not corresponding to an extraso¬ 
lar planet. 

2.2. HD1461 

Diaz et al. 2015 (in prep.) analyzed the HARPS RV 
data of HD 1461 and found the presence of two plan¬ 
ets orbiting at 5.77 and 13.50 days with 0.05 and 0.09 
of eccentricity, respectively. In addition, the star also 
shows a long-term magnetic cycle in the variation of 
its calcium activity index. The effect of the magnetic 
cycle is seen in the RVs of HD 1461, and Diaz et al. 
2015 (in prep.) fit it using a Keplerian with a period of 
3508 days and an eccentricity of 0.10. 

After fitting a three Keplerian model to the RVs of 
HD 1461 with parameters fixed to these previous val¬ 
ues, the RV residuals show a signal at a one-year pe¬ 
riod (Figuretop right panel). This signal has an am¬ 
plitude of 1.31 ± 0.20 ms“^ when a model with three 
Keplerians plus a one year sinusoid is fitted to the data. 

As for HD 128621, the anti-correlation between the 


^ one pixel on the HARPS CCD represent 820 m s ^ 
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Fig. 1.— Generalized Lomb Scargle (GLS) peri- 
odograms ( [Zechmeister & Kurster|2009) l of the RV 

residuals of our four stars after removing their plane¬ 
tary signals and their magnetic cycle effect (see text in 
Sec 1^ for more details about the fitted models). The 
red horizontal lines represent a false alarm probabil¬ 
ity level of 10%. For each star the one-year signal is 
significant and is highlighted by the red regions. 
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Fig. 2.— Top: RV residuals of HD128621 in blue af¬ 
ter fitting a 3rd order polynomial. The red curve 
corresponds to the fitted one-year signal. Bot¬ 
tom: B ary centric Earth radial velocity (BERV) of 
HD 128621, which is the velocity of the Earth in the 
direction of HD 128621. 


BERV for HD 1461 and the one year signal seen in the 
RVs of the star shows that the two signals are in op¬ 
position of phase (Figuretop right panel). Therefore 
the revolution of Earth around the Sun seems to be de¬ 
tected in the RVs of HD 1461. 

2.3. HD154088 


Mayor et~ST] ( |2011| ) announced the detection of a 
super-Earth orbiting around HD 154088. After fitting 
a Keplerian model with the planet parameters reported 
in |Mayor et al. ( 2011| l and a 2nd order polynomial drift 
to account for the magnetic cycle, the RV residuals 
show a significant signal at a one-year period (Figure 
bottom left panel). This signal has an amplitude of 
0.76 + 0.16 m s“^ when a model with one Keplerian, a 
2nd order polynomial, plus a one year sinusoid is fitted 
to the data. 

As for HD 128621 and HD 1461, the anti-correlation 
between the BERV for HD154088 and the one year 
signal seen in the RVs suggests that the signal of 
Earth orbiting the Sun is detected (Figure [^bottom left 
panel). 

2.4. HD31527 

A planetary system of three planets orbiting HD31527 
has been announced by [Mayor et al.| ( [2()ll| ). The ac¬ 
tivity index variation for the star is very low, with an 
average logiR' hk) level of -4.96. We therefore do not 
expect any impact of activity on the RVs. When fit¬ 
ting the three planets with the parameters reported in 
Mayor et al.| ( |2Q11| ), a signal with a period close to 
one year is seen in the RV residuals (Figure bottom 
right panel). This signal has an amplitude of 0.59 + 
0.15 ms“^ when a model with three Keplerians plus a 
one year sinusoid is fitted to the data. 

As for the other stars studied in this paper, the op¬ 
position in phase between the one-year signal detected 
in the RV residuals and the BERV for HD31527 sug¬ 
gests that the signal of Earth is detected in the stellar 
RVs (Figure [^bottom right panel). 

3. Origin of the one-year signal 

We show in the preceding section that HARPS RV 
data of several stars that have been intensively ob¬ 
served during several years are affected by a one-year 
signal that is in opposition of phase with the yearly 
variation of the BERV of the target. The same oppo¬ 
sition in phase observed in the RV of several stars 
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Fig. 3.— Correlations between the RV induced by the one year signal and the barycentric Earth radial velocity (BERV) 
for our four stars in red. The blue points represent the underlying data with error bars. 


cannot be a coincidence. It suggests that the signal 
that we are detecting is due to the Earth orbiting 
the Sun. Several different effects could be responsi¬ 
ble for this annual signal: 


imprecision in the code to calculate the BERV 


([Wright & Eastman|2014^ [Kopeikin & Qzer- 

|noy|1999> , 


• incorrect stellar coordinates, 


• telluric absorption line contamination, 

• incorrect observatory coordinates, 

• incorrect reduction of observing times, 

• uncorrected systematics in seasonal varia¬ 
tions at the telescope. 


and others. It is difficult to rule out an effect by 
studying a group of stars, because each star will 
have a different BERV. However, this is not the goal 
here and we will demonstrate that the yearly signal 
is, in the presented example, mostly created by im¬ 
perfections of the CCD that induce a deformation 
of spectral lines when the stellar spectrum is red- 
and blue-shifted on the CCD due to Earth orbiting 
the Sun. 


3.1. Measuring the radial velocity of individual 
spectral lines 

When Earth is revolving around the Sun, the BERV 
of a target is changing as a function of time, and this 
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Fig. 4.— RV residuals as a function of time for 
HD 128621 derived from the HARPS spectral order 61, 
after fitting a 3rd order polynomial to the data. The 
one-year signal is significant in the corresponding GLS 
periodogram (red region). 


translates into a shift of the recorded spectrum on the 
CCD due to the Doppler effect. If a spectral line passes 
over an imperfection of the CCD, this spectral line can 
undergo a deformation that will introduce a RV varia¬ 
tion correlated or anti-correlated with the BERV. 

To test this hypothesis, we first looked at the RV 
of each spectral order individually, as some spectral 
orders might be more affected than others. In Figure 
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we can see the RVs derived using only the spectral 
order 61. This spectral order spans a wavelength be¬ 
tween 6145 and 6215 A. Despite the absence of rel¬ 
evant telluric lines, the RV of this spectral order is 
nevertheless strongly affected by a one-year signal. 

To go a step further and detect which lines in spec¬ 
tral order 61 are responsible for the yearly signal, we 
fitted each individual line stronger than 20% in depth 
with a Gaussian and selected the center of the Gaussian 
as the average wavelength of the line. This procedure 
is performed on all the stellar spectra to get the aver¬ 
age wavelength position of each line as a function of 
time. The variation in average wavelength position as a 
function of time is then translated into a RV variation. 
Figure shows the RV variation of four spectral lines 
as a function of time. The top two lines are very close 
in wavelength, however one of them shows a strong 
one-year signal, while the other does not. The same 
behavior is observed for the two spectral lines at the 
bottom of Figure The two lines not affected by 
the yearly signal show no correlation between RV and 
BERV (see Figure. [^, while the two other lines show a 
strong anti-correlation. This anti-correlation is similar 
to the one observed when using the RV derived with 
all the spectral lines in the spectrum (see Figure [^, 
and therefore those two lines are examples of lines 
contributing to the detected yearly signal. 

3.2. Origin of the CCD imperfections deforming 
spectral lines 

The HARPS CCD is composed of a mosaic of 
two 2x4K EEV44-82 CCDs (HARPS detector final 
design revie\\0). Each CCD device is subdivided 
in 1024x512 blocks of pixels (1024 in the cross¬ 
dispersion direction and 512 in the spectral direction) 
with a pixel size of 15 yum within each block. The 
block boundaries show sometimes a 1 % quantum ef¬ 
ficiency (QE) variation over 1 row or column due to 
photolithography stepper mismatches. This variation 
in QE is induced by a variation in pixel size at the area 
boundaries, due to difficulties to perform the stitching 
of all the blocks together. These boundaries will be 
called block stitchings hereafter. 

To get a wavelength solution for each pixel illumi¬ 
nated by both spectrograph fibers, a Thorium-Argon 
lamp is used to illuminate each fiber. The emission 
lines in the Thorium-Argon lamp create a finger print 


^ http ://www. eso. org/sci/php/optdet/instmments/harps/FDR-Harps- 
detectors.pdf 


on the CCD, and by knowing precisely the wavelength 
of each individual line, it is possible to derive a wave¬ 
length for all the pixels in each spectral order. Because 
Thorium-Argon emission lines are sparse on the de¬ 
tector, a polynomial fit is performed on each order 
to get a wavelength for all pixels. The wavelength 
solution is correct if all pixels have the same size, 
which is not the case for the boundaries between the 
1024x512 blocks constituting the CCD. This poly¬ 
nomial fit therefore introduces an error in wave¬ 
length solution at those boundaries. This effect has 
been shown by Wilken et al. ( 2010| ) (see Eigure 4 of 
their paper), where the authors compare wavelength 
solutions obtained with a laser frequency comb and a 
Thorium-Argon lamp. Every 512 pixels in the spec¬ 
tral direction, a step as high as 40 ms“^ (l/20th of 
a pixel on HARPS) is observed when comparing the 
two wavelength solutions. These differences can be 
explained by the larger number of regular emission 
peaks present in the laser frequency comb that can re¬ 
solve the difference in pixel size at the block stitch¬ 
ings, which the spatially sparse Thorium-Argon emis¬ 
sion lines can not do. 


In Eigure [ 7 ] we show the part of the spectrum 
around both sets of spectral lines studied in the pre¬ 
ceding section (6154.68 and 6155.56 A, and 6198.74 
and 6199.86 A). The two different spectra in blue and 
red correspond to observations done at the extremes of 
the BERV, highlighting the shift of the spectrum on the 
CCD. As we can see, the spectral lines at 6155.56 and 
6199.86 A affected by the yearly signal cross the block 
stitchings represented by the green vertical dashed line 
when Earth orbits the Sun. The two other lines do not 
cross those boundaries and therefore do not show the 
one-year variation. 


4. Correcting for the one-year signal 

In the preceding section, we demonstrated that 
the boundaries between the blocks constituting the 
CCD, so-called block stitchings, are responsible for 
the yearly signal detected in HD128621, HD1461, 
HD154088 and HD31527. A simple solution to get 
rid of this yearly signal is to remove the spectral lines 
passing over these boundaries from the correlation 
mask used to derive the RVs. Because each star has 
its own gamma velocity, the position of the spectral 
lines on the CCD will differ from one star to the other, 
and therefore a different correlation mask has to be 
generated for each star. 
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Fig. 5.— RV residuals as a function of time for HD 128621 derived from different spectral lines present in the HARPS 
spectral order 61, after fitting a third order polynomial to the RV of each line. The corresponding GLS periodograms 
are also shown and highlight in red the presence of the one-year signal for spectral lines at 6155.56 and 6199.86 
A, and the absence of the signal for the nearby lines at 6154.68 and 6198.74 A. 
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Fig. 6.— Correlations between the RV derived from different spectral lines present in the HARPS spectral order 61 and 
the BERV for HD 128621. The red lines represent a linear fit to the data, and the corresponding correlation coefficient 
is reported in the legend. A strong anti-correlation is observed for spectral lines at 6155.56 and 6199.86 A 
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Fig. 7.— Zoom on small spectral regions of spectral order 61. The blue spectrum corresponds to an observation of 
HD 128621 when the BERV was maximum and the red spectrum when the BERV was minimum. The green dashed 
vertical lines corresponds to the position of the block stitchings. Spectral lines at 6155.56 and 6199.86 A showing the 
yearly effect cross these block stitchings, while the other lines, not affected, do not. 
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To select the spectral lines at wavelength Ai that 
cross the block stitching when Earth is revolving 
around the Sun, the first step is to shift their wave¬ 
length to the stellar gamma velocity y*. The new 
wavelength in the stellar rest frame T*,/ can be com¬ 
puted using the Doppler effect formalism = 

Ai ^ - ij. The second step is to identify the 

wavelength Absj of the block stitchings, which can 
be done using the wavelength solution and keeping in 
mind that the block stitchings appear every 512 pixels 
in the spectral direction. Finally, all the spectral lines 
in the correlation mask that have a wavelength in the 
stellar rest frame obeying to the relation 


'^BSj - i^inBERV -W< Ai^j < Absj + maxBERV + fk (1) 

should be excluded. Here W corresponds to the width 
of a spectral line. We fixed this value to 10kms“^ to 
exclude the wings of each spectral line[^ 

In Figure we illustrate the process of excluding 
the spectral lines from the correlation mask taking the 
example of spectral order 61. The top panel shows 
two spectra of HD 128621 taken at the extremes of the 
BERV. Each block stitching is localized in the wave¬ 
length domain, and the blue region across each block 
stitching correspond to the region where spectral lines 
should be excluded. The bottom panel shows a spec¬ 
trum of HD 128621 taken when the BERV was close 
to zero. All the spectral lines appearing in the correla¬ 
tion mask and shifted to the stellar gamma velocity are 
shown, and the ones falling close to the block stitching 
are excluded (red lines). Note that on the blue side of 
the spectral order, a line is exclude while not falling 
on a block stitching. There is some overlap between 
the spectral orders, and this line is rejected in spectral 
order 60. 

For each star, two different correlation masks are 
generated. One with only the lines falling close to the 
block stitchings, and another one without those lines. 
The RVs measured with both correlation masks can 
be seen in Figure It is clear that the spectral lines 
crossing the block stitching are the origin of the one- 
year signal (see left panel). Once those spectral lines 
are removed from the correlation mask, the yearly sig¬ 
nal disappears (see right panel). Note that this Figure 
should be compared with Figure 

^The average full width at half maximum of a line is a few km for 

a solar like star rotating slowly jGray|200^ 


5. Conclusion 


We have detected a one year signal in the data of 
several stars intensively observed with HARPS dur¬ 
ing several years. This signal is induced by spectral 
lines crossing the block stitchings of the HARPS CCD 
when the stellar spectrum is shifted back and forth on 
the detector due to the Earth’s orbit around the Sun. 
At the position of each block stitching, the wavelength 
solution is not correct because the sparse emission 
spectrum of the Thorium-Argon lamp cannot resolve 
the different pixel size. Spectral lines crossing block 
stitchings therefore undergo a deformation when pixel 
position is translated into wavelength using the wave¬ 
length solution, which induces a RV variation. Remov¬ 
ing those spectral lines from the correlation mask used 
to derive the RV with the cross-correlation technique 
allows to get rid of this yearly signal. Note that a differ¬ 
ent correlation mask for each star has to be generated 
due to a different stellar gamma velocity, and therefore 
a different position of the spectral lines on the HARPS 
CCD. 


The amplitude of the one-year signal is different 
for each star, because of different BERV and different 
spectral lines crossing the block stitchings of the CCD 
(see Figure [^. If the amplitude of the BERV is only 
a few km s” (maximum is 30kms“^), then the stel¬ 
lar spectra will only be shifted on the CCD by a few 
pixels, which will not deform enough the spectral lines 
crossing the block stitchings and therefore not induce a 
significant signal at one year. In addition, if only weak 
spectral lines, that only have a small weight in the fi¬ 
nal RV computation as described in |Pepe et al.| ( |2003| , 
cross block stitchings, the induced RV variation will 
also be negligible. The combination of these two ef¬ 
fects explains why some stars are strongly affected 
by the one-year signal, while others are effected to 
a lesser degree. 


The number of spectral lines rejected in the correla¬ 
tion mask is of the order of 20%. This can therefore be 
a problem when observing faint targets, for which the 
signal-to-noise-ratio of the spectra is low. A possible 
way of removing the yearly signal without rejecting 
spectral lines would be to include as free parameters 
the pixel sizes at all block stitching, and thus derive a 
corrected wavelength solution (Coffinet et al. 2015, in 
prep.). 

For calibration sources such as Fabry-Perot 


etalons or laser frequency-combs (Eangellier et al. 


2014[|Wilken et al.|2012| ) that have a higher density of 
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Fig. 8.— Top: The blue spectrum corresponds to spectral order 61 of an observation of HD 128621 when the BERV 
was maximum and the red spectrum, spectral order 61 when the BERV was minimum. The green dashed vertical 
lines correspond to the position of the block stitchings and the blue regions to the zone where spectral lines should be 
excluded. Bottom: Spectrum corresponding to order 61 of an observation of HD 128621 when the BERV was close 
to zero. The dashed vertical lines are the spectral lines present in the correlation mask, shifted to the stellar gamma 
velocity. The spectral line falling inside the excluding regions are highlighted in red and should be excluded from the 
correlation mask when deriving RVs to prevent the one-year block stitching effects. 
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Eig. 9.— Left: GLS periodograms of the RV residuals of our four stars after removing their planetary signals and 
their magnetic cycle effect (see text in Sec for more details about the fitted models). The RVs are derived using 
the correlation masks with only the spectral lines crossing the block stitchings. The periods close to one year are 
highlighted by the red regions, and are significantly above the 10% false alarm probability level (red horizontal lines). 
Right: Same plots but for the correlation masks cleaned from spectral lines crossing the block stitchings. 
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emission lines on the detector, it is possible to detect 
the CCD block stitchings and therefore directly obtain 
a correct wavelength solution. In addition to a better 
RV precision, these calibration sources should pro¬ 
duce data that are not affected by the one-year signal 
described in this paper. 

This yearly signal is extremely problematic when 
looking for long-period small-amplitude planets as it 
will inject signal at 365.25 days and the first cor¬ 
responding harmonics, generally 182.63 and 121.75 
days. It is therefore mandatory to correct this yearly 
signal from the HARPS RV data if such planet want to 
be detected. For a user detecting a spurious signal 
at one year in HARPS RV data, two different meth¬ 
ods can be used to suppress this signal. Using the 
pipeline, it is possible to recalculate a new correla¬ 
tion mask based on the technique described in Sec¬ 
tion 4 and reprocess the CCF routine. Otherwise 
fitting a sinusoidal signal with a one-year period is 
precise enough to a few dozens of cms“^. A care¬ 
ful check of the anti-correlation between the fitted 
signal and the BERV should however be done, to 
prevent removing one-year planetary signals that 
could exist in the data. 

We thanks the referee for his valuable comments 
and suggestions that improved the first version of 
the paper. X.D. thanks NASA’s Transiting Exo¬ 
planet Survey Satellite (TESS) mission for partial 
support via subaward 5710003554 from MIT to 
S AO. We are grateful to all technical and scientific col¬ 
laborators of the HARPS Consortium, ESO Headquar¬ 
ters and ESO La Silla who have contributed with their 
extraordinary passion and valuable work to the success 
of the HARPS project. 
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